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AWARDS ABSTRACT 

The invention relates to providing periodic gain adjustment in plants 
of irreducible order, n, as a way of modulating control energy to higher 
frequencies and making the harmonics add in such a way that a nonminimum- 
phase system appears to be minimum phase using a technique based on the 
notion of a mathematical "lifting" of extended horizons in which a se- 
rial-to-parallel conversion is performed on N>n liftings taken from the 
control input and output signals, and mappings are considered between the 
vectorized quantities. This lifting technique is also applicable to com- 
munications for equalization of a nonminimum phase channel. 

Referring to FIG. 1, an improved method for providing periodic gain 
adjustment in plants of irreducible order, such as for control of 
electrical, mechanical or chemical systems or for equalization of 
communication channels, comprises the choosing of a horizon time N>n of 
liftings in periodic input and output signal windows and py, respec- 
tively, where N is an integer chosen to define the extent (length) of each 
of the windows p„ and Py, and n is the order of the irreducible input/output 
plant model given by an equation of the form 
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where polynomials A and B are relatively prime, such that the horizon time 
for liftings have a total length N=m+£+p+q+2n-l , and in which the number 
of elements in the liftings taken from the input signal is greater than 
the number of elements Cy in the liftings taken from the output signal for 
output tracking in control systems and vice versa (a„<a ) for input tracking 
in communication systems. The extended horizon liftings for output 
tracking over a selected input signal window p„ and a selected output 
signal window py are of the form 
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where jn, 2, p, n, q, and n-1 are elements of the lifting, and m^O and 2^0 
are arbitrary, q=0, p e is an arbitrary (or null) 0-X vector chosen 
identically in both and Py and n>0 is the order of the irreducible plant. 
If the system is obtained Xjy zero-order hold digitization of a continuous- 
time plant, the integer q^O can be chosen arbitrarily. Since there are m+2 
zeros on the left in Py and only m zeros on the left in p„, the number of 
elements is greater than Oy for output tracking. For input tracking in 
a communications system, the extended horizon liftings are of the same 
general form, but now there are m+2 zeros on the left in p„, and only m 
zeros on the left in Py so that the number of elements is less than Cy. 

FIG. 2 is a block diagram of a generalized lifting system model 
Gfpu,py) applicable to control systems in general. FIG. 3 illustrates 
graphically the determination of a small matrix from input and output 
signal liftings in respective windows p„ and py. FIG. 4 is a block diagram 
of a generalized lifting system model under zero annihilation conditions BSu =0 
and . FIG. 5 is a block diagram of a zero annihilation periodic 
(ZAP) control law having a closed-loop system arising from the output 
tracking of extended horizon liftings. FIG. 6 is a pole-zero plot of a 12- 
state nonminimum phase plant with a sampling interval T=0.025sec in which 
poles are represented by "X" and zeros by "0." FIG. 7 is a graph of open- 
loop output response of the plant in FIG. 6 to an initial condition. FIGs. 
8a and 8b are graphs of respective output and input closed-loop response 
of the plant in FIG. 6 using deadbeat control with Lozano's 2n-lifting 
(a„=Oy) with a vibration damping objective. FIGs. 9a and 9b are graphs of 
respective output and input response using ZAP control of FIG. 5 with an 
extended horizon input lifting defined by m=0, 2=40, p=0 , n=12 , q=0 , 
N=63. In this case a„=52 and Oy=X2, so that a„>ay. FIG. 10 illustrates 
application of the present invention to communications for equalizing a 
nonminimum phase channel. 

The novelty of the invention resides in the use of extended horizon 
(window length) for periodic gain adjustments in control systems with 
output tracking and for equalization of communication channels with input 
tracking, coupled with other conditions that are similar for both output 
tracking control and input tracking equalization of communication channels. 
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EXTENDED HORIZON LIFTINGS FOR PERIODIC GAIN ADJUSTMENTS IN 
CONTROL SYSTEMS, AND FOR EQUALIZATION OF COMMUNICATION CHANNELS 


ORIGIN OF INVENTION 

The invention described herein was made in the performance 
of work under a NASA contract, and is subject to the provisions 
of Public Law 96-517 (35 USC 202) in which the contractor has 
5 elected not to retain title. 


TECHNICAL FIELD 

The invention relates to providing periodic gain adjustment in 
control design as a way of modulating control energy to higher fre- 
10 quencies and making the harmonics add in such a way that a nonmini- 
mum— phase system appears to be minimum phase using a technique 
based on the notion of a mathematical "lifting" in which a se- 
rial-to-parallel conversion is performed on the control input and 
output signals, and mappings are considered between the vectorized 
15 quantities. This lifting technique is also applicable in communi- 

cations for equalization of a nonminimum phase channel. 

BACKGROUND OF THE INVENTION 

In the prior art, nonminimum-phase systems could only be 
20 controlled by using low bandwidth, sluggish controllers. In lay- 

man's terms, nonminimum-phase systems move in a direction initially 
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opposite to the direction they are pushed. Such systems arise 
when, for example, one attempts to dampen vibrations on a flexible 
structure in optical instruments by acting at a location remote 
from the flexible structure in the case of a "noncolocated" control 
system. 

Nonminimum-phase systems arise in many aerospace applications 
of noncolocated control, such as spacecraft slewing, instrument 
pointing/ tracking, flexible robotics, acoustic systems, chemical 
process control systems, compact disk controllers, floppy disk 
controllers and even real-time equalization of nonminimum-phase 
communication systems. An object of this invention is to achieve 
high performance control for the class of nonminimum-phase systems 
that encompasses these and other applications which are referred to 
generically as "plants." 

A Zero Annihilation Periodic (ZAP) control law was introduced 
by the present inventor, David S. Bayard, in publications titled 
"Globally stable adaptive periodic control," Jet Propulsion 
Laboratory, Internal Document JPL D-9448, February 3, 1992 and 
"Zero annihilation methods for direct adaptive control of 
nonminimum-phase systems," Proc. Seventh Yale Workshop on Adaptive 
and Learning Systems, Yale University, May 1992, for controlling 
nonminimum-phase systems using stable inversion of the plant. The 
general ZAP approach is based on the notion of a mathematical 
"lifting" performed on the plant input and output signals with 
mapping between the vectorized input and output quantities. 
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A key property of the ZAP liftings which makes them so useful 
is that the transmission zeros of the lifted and vectorized plant 
input signal are annihilated (i.e., placed to the origin). This 
ZGro~annihilation (ZA) property allows the vectorized plant signals 
to be stably inverted using standard control methods. The result 
is important to many areas of control, communications, and signal 
processing where a stable plant inverse is often desired but not 
possible due to nonminimum-phase restrictions. 

In a paper by R. Lozano-Leal, "Robust adaptive regulation 
without persistent excitation," IEEE Trans. Automatic Control, Vol . 
34, pp. 1260-1267, December 1989, a multirate sampling method is 
presented which allows stable inversion of any linear time-invar- 
iant finite-order plant signal. Not surprisingly Lozano's lifting 
has also been applied to developing stable adaptive control algo- 
rithms for nonminimum phase systems. [Bayard, supra and Lozano- 
Leal supra']. See also R. Lozano et al., "Singularly-free adaptive 
pole placement using periodic controllers," IEEE Trans. Automatic 
Control, Vol. 38, pp. 104-108, January 1993. 

in this invention is an extension of the prior art 
represented by Bayard, supra, and Lozano-Leal, supra, to "extended 
horizons," i.e., to extended periodic windows over which liftings 
(samples) are taken in adaptive controllers of the type described 
by Lozano et al., supra, for example. Extended horizon liftings 
are crucial for control gain reduction in order to allow practical 
implementations of control in nonminimum— phase systems. The dis— 
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tinction between extended horizon liftings of the present invention 
and the prior art will now be described. 

Since Lozano~Leal ' s liftings are distinct from the present 
invention in that it utilizes a horizon size of iV=2n, where n is 
5 the plant order and N is the length of the window for the liftings, 
and the number of elements in the system input equals those in the 
system output (i.e., a„=Oy) , it will be denoted as the "2n-lifting, " 
while in the present invention N>n and preferably N>2n-1 in a large 
class of systems in which o^>Oy for output tracking (mapping) , such 
10 as in feedback control systems, and cr„<ay for input tracking, such 

as in communications systems, and which enjoy the same zero anni- 
hilation properties as N=2n. Unlike the special case of Lozano- 
Leal 's and other prior art liftings, the present invention allows 
the use of all extended horizon liftings with N greater than n and 
15 preferably equal to or greater than 2n-l with CT„>ay for output 
tracking (OT) and o^<Oy for input tracking (IT). An important 
consequence is that such extended horizon liftings lead to plant- 
-inverse controllers with significantly reduced control gains. 
This overcomes a problem associated with Lozano-Leal 's lifting, 
20 where N=2n and a^=ay, which has prevented its use in many appli- 
cations of practical interest that the present invention will 
reach. To illustrate the present invention, a simulation example 
is provided below in which the peak control requirejnent is rcducod 
by four orders of magnitude using extended horizon liftings. 
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It will also be shown that as a dual result, a related class 
of extended horizon liftings enables equalization of nonminimum- 
phase channels in communication systems. This overcomes the stan- 
dard problem of inverting the channel in a stable fashion. In this 
xnvention, the extended horizon property allows channel inversion 

by least squares estimation, which provides smoothing in the case 
of noise. 

As background information, consider a plant input/output 
model , 


A(z-i)=l+£a,z-; 

-1=1 i=l 

where the polynomials A and B are assumed to be relatively prime. 
It is assumed that , so that the polynomial B can be factored 
uniquely into the form B(s ') = 2 -‘'b^fl(s ■) where fl (2 ■) is a monic 
polynomial and d=l is the plant delay. The choice d=l is lor 
convenience only and is not a fundamental restriction. In the case 

that dal , all subsequent expressions can be appropriately modified 
without loss of generality. 
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Choose some horizon time N>n. The system of Eq. (1) is 
iterated to give the following system of linear equations, 

Y(k)=A^Y(k) +A2Y{k-l)+B^U{k) +B2U(k-l) (2) 
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where , 
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( 3 ) 


Y(k) = plant output 
15 U(k) = plant input 

Aj = lower triangular Toeplitz, with first column 

• • • , Of m • .,0]^ 

Aj = upper triangular Toeplitz, with first row 

[Of ~3j 

20 J3i = lower triangular Toeplitz, with first column 

[ f • • • r fOf . . . 0 ]^ 

JBj = upper triangular Toeplitz, with first row 
[ 0 , • • • fO fb„f m m m fb2^ 


25 
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Example 1 Let n=3 and iV=4 . Then Eq. (2) becomes. 


YAk-3 

YAk-2 

74*^1 [0 -Oe -«2 -«ll 0 0 0 

Yik^2 0 0 -03 -02 1 -«! 0 0 0 

74^.3 % 0 0 - 03|-«2 -«1 0 0 

74 ;, J [o 0 0 0 |-a 3 -«2 -«1 oj 

V4k^2 

Ytk*3 

Yik*A 


0 0 Jb3 bj I bi 0 0 0 

15 0 0 0 b3 I bj bj 0 0 

^ 0 0 0 0 I b3 bj bi 0 

0 0 0 0 I 0 b3 b2 bi 


20 It is convenient to combine terms involving Y(k) in Eq. (2) and 
rearrange to give the following lifting of P. Albertos, "Block 
multirate input-output model for sampled-data control systems," 
IEEE Trans. Automatic Control, Vol. 35, No. 9, pp. 1085-1088, 
September 1990. 



25 
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Albertos' Lifting: 

Y{k) =AY{k-l)*HU{k) +BU{k-l) 


(4) 


where , 


5 


A = (I-A^)~^A2 


(5a) 



(5b) 

B = (I-A^)% 

(5c) 


10 

It is noted that since Aj is lower triangular with zeros on the 
diagonal, the quantity (I-Aj) is always invertible. Hence the quan- 
tities in Eq. (5) always exist. 

Polynomial A is divided into B to give impulse response 

15 sequence {hi} , 


A{z-^) 


00 



( 6 ) 


The Markov parameter sequence {h^} is not assumed to be convergent 
20 (i.e., the system may be unstable). Using the Toeplitz structure of 

and Bi and relation, (Eq. (6)), it can be shown that the matrix 
H in Eqs. (4) and (5) can be written in terms of the impulse 
response parameters , 
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( 7 ) 


This is the desired expression for H, i.e., H = lower triangular 
Toeplitz, with first column since the delay is 
unity by assumption (i.e., d=l) , the matrix H has a nonzero diag- 
onal (i.e., hi?fO), and is always invertible. 


STATEMENT OF THE INVEHTTOM 

In accordance with the present invention, an improved method 
for providing periodic gain adjustment in plants of irreducible 
order, such as for control of electrical, mechanical or chemical 
systems or for equalization of communication channels, comprises 
the choosing of a horizon time N>n of liftings in periodic input 
and output signal windows p. and p„ respectively, where ff is an 
integer chosen to define the extent (length) of each of the windows 

P. and p„ and n is the order of the irreducible input/output plant 
model given by an equation of the form 


A(z-i)=l+^a^z-^; B(z-i)=|^b.2-i 

- 1=1 


where polynomials A and B 


are relatively prime, such that the 
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horizon time for liftings have a total length fptgt2n-l , and 

in which the nuntoer of elements o. in the liftings taken from the 
input signal is greater than the number of elements o, m the 
liftings taken from the output signal for output tracking in 
control systems and vice versa (o.<o,) lor input tracking in 
communication systems. The extended horizon liftings for output 
tracking over a selected input signal window n. and a selected 
output signal window py are of the form 





m 


p =[^.., 0 , o,...,o, p b,...,o, i,i,....il 



where m, f, P, n, g, and n-1 are elements of the lifting, and mad 
and «0 are arbitrary, g=0, p s is an arbitrary (or null) 0-1 
vector chosen identically in both p, and p, and n>0 is the order of 
the irreducible plant. If the system is obtained by zero-order 
hold digitization of a continuous-time plant, the integer gaO can 
be chosen arbitrarily, since there are m+f zeros on the left of p, 
and only m zeros on the left of p., the number of elements c. is 
greater than c, for output tracking. For input tracking, the 
extended horizon liftings are of the same general form, but now 
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there are m+£ zeros on the left of p„, and only in zeros on the left 
of Py so that the number of elements is less than Oy. 

5 BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 graphically depicts partial horizon vectors Ug and ¥g 
using windows p„ and Py. 

FIG. 2 is a block diagreun of a generalized lifting system 
model G(pg,py) defined by Equation 19 applicable to control systems 
10 in general. 

FIG. 3 illustrates graphically the determination of a small 
matrix from input and output signal liftings in respective win- 
dows Pu and Py. 

FIG. 4 is a block diagram of a generalized lifting system 
15 model under zero annihilation conditions BSj=0 and A(Sy)^=0 defined 
by Eqs . (27) and (28). 

FIG. 5 is a block diagreun of a zero annihilation periodic 
(ZAP) control law having a closed-loop system arising from the out- 
put tracking of extended horizon liftings of Equations (35) and 
20 (36). 

FIG. 6 is a pole-zero plot of a 12-state nonminimum phase 
plant with a sampling interval T=0.025sec in which poles are repre- 
sented by "X" and zeros by "0." 

FIG. 7 is a graph of open-loop output response of the plant in 
25 FIG. 6 to an initial condition. 
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FIGs. 8a and 8b are graphs of respective output and input 
closed-loop response of the plant in FIG. 6 using deadbeat control 
with Lozano's 2n-lifting {a^=^Oy) with a vibration damping objective. 
FIGs. 9a and 9b are graphs of respective output and input 
5 response using ZAP control of FIG. 5 with an extended horizon input 

lifting defined by m=0, 2=40, p=0, n=12 , q=0 , N=63 . In this case 
0^=52 and Oy=12, so that a„>ay. 

FIG. 10 illustrates application of the present invention to 
communications for equalizing a nonminimum phase channel. 

10 

DETAILED DESCRIPTION OF THE INVENTION 

A class of liftings {N>n) will first be defined by generaliz- 
ing the lifting of Albertos (Eq. (4)). For this purpose, it will 
be useful to construct the "small" vector Y„(k) from Y(k) as 
15 follows, 

YJk)&S^Y(k)eR°y 

where is a selection matrix which sifts out Oy elements of 

Y(k) for inclusion into Yg(k ) . 

The matrix Sy is most conveniently constructed from a 0-1 vec- 
20 tor py. For example, if py = [1, 1,0,0] then Y^(k) e contains the 
first two elements of Y(k) e R* and the selection matrix is given 


f 


10 0 0 
0 10 0 


as 
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A systematic method to construct S, is defined as follows: form a 

diagonal matrix from the entries of rector and then remove all 
rows made up entirely of Os. This construction defines the mapping 
W : for which one can write s, = W(p,) . since p, can be 

5 uniquely reconstructed by a logical "or" over the columns of S„ 
the mapping W is one-to-one. 

Using the above notation, the following "small" vectors are 
defined. 


10 

YJk)&SyY(k); 


(8a) 


UAk)&S^U(k); 


(8b) 


Yf{k)^S;Y(k)j 

Sy^{p^)eR 

(8c) 


15 


where p, and p. are specified 0-1 window vectors, and is defined 
as the 0-1 complement of p,. The vector rf(Jc) in Eq. (8c) is de- 
noted as the complementary output since it comprises all elements 
of the vector Y(k) which are not included in ¥,(k) . 

A formula to reconstruct r(k) from r.(k) and y/(*) is now 
derived. First, it is to be noted that the quantity ^^-[Sy , 

is a permutation matrix. Hence its Inverse is its transpose, i.e., 
r which upon expanding gives 


25 




( 9 ) 
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Multiplying both sides of Eq. (9) on the right by Y(k) gives the 
desired formula. 




( 10 ) 


Example 2. FIG. 1 graphically depicts partial horizon vectors 
y. and r, for the limited liftings case of n=3 and N=6 , as deter- 
mined by window vectors p,=[0, 1 , 1, 1, 0, 0] and p,=[0,0, 1, 1, 1,0] . In 
this case, py‘"=[l,l, 0, 0,0,1] and one can compute. 


Su=^{Pu)= 


0 1 0 0 0 0 

0 0 1 0 0 0 ; 

0 0 0 1 0 0 


Sy=W{p^) 


0 0 1 0 0 0 
0 0 0 1 0 0 
0 0 0 0 1 0 




1 0 0 0 0 0 
0 1 0 0 0 0 
0 0 0 0 0 1 


In the present invention, a new family of extended liftings 
{N^2n-1) will be defined using the notation developed above for 
Albertos' «-lifting in Eg. (4). As a k^y step. It will be assumed 
that U(k) is chosen as zero outside the window defined by 
Mathematically this can be written as. 


{I-sJsjU(k)=Q ( 11 ) 

Consider the following nonminimal state-space realization of 
Albertos' iV-lifting Eg. (4) determined by 
in Eqs. (8) and identity Eq. (11), 


using the small vectors 
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Y(k)- 

A BS^ 

Y{k-l) 

\hsJ 

Usik) ' 

0 0 

Usik-1) 

h 

j 


Us(k) 


( 12 ) 


5 







(13) 


Substituting for the state in the output Eq. (13) gives the alter- 
native output equation, 
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Y,(k)=[SyA 


SyBSj] 


Y(k-l) 


^SyHSjujk) 


(14) 


15 


For analysis purposes, it is convenient to transform the 
open-loop plant using the similarity transformation T where. 


Ysik) 

Ys'ik) 

U,{^) 


Y{k) 

Usik) 


(15) 




C 


Sy 10 

20 


M 

; C = [Sy 0]; M- 

y 1 

0 |J 


S;^W{P^y) (17) 

It can be verified that transformation T is square and invertible. 
25 Furthermore, T is in the form of a permutation matrix which 
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reorders the state such that the components of Y, appear first, and 
the remaining elements follow in the specified order. Since T is a 
permutation matrix, the inverse of T is given simply by its trans- 
pose (cf., Barnett Afa trices; Methods and Applications, Clarendon 
Press, Oxford, England, 1990, pp. 374), i.e., 

( 18 ) 

Transforming the open-loop dynamics of Eq. (12) by the similarity 
transformation of Eqs. (15) and (16), gives rise to a very useful 
representation denoted as the Generalized Lifting System Model, 


Generalized Lifting System Model, G(p^,p ): 




Ysik) 


Usik) 



SyAS^ Sy4(s/)^ 

SyASy SyA{Sy)'^ S/BS^ i " 1 ) 

0 0 0 


SyHSj 

SyHSj 


Us(k) 


( 19 ) 


It is noted that the generalized lifting G(p„, p^) of horizon 
length N<n is defined uniquely by the choice of selection windows 
p„ and py, from which the matrices are calculated. The 
generalized lifting system model of Eg. (19) is depicted in the 
block diagram of FIG. 2. It is seen that r. and r/ form two 
coupled subsystems which are driven by a common input U^. 


It is 
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also noted that the transmission zeros of the transfer function 
from Ug(k) to Yg(k) are affected by the choice of windows Py* 
This is a key feature which will be used to advantage in 
applications of the present invention. 

The lifting G(p„, Py) generalizes a number of existing results. 
For example, Albertos' lifting in Eq. (4), taken from the paper of 
Albertos, supra, is equivalent to a choice of periodic windows 
(horizon length N) as shown in Eq. (20) for input and output 
liftings without regard to the plant order which should a priori be 
excluded for consideration in applications of the present 
invention . 


N 


( 20 ) 




N 


Py~[lflf ■ • • 


\ 

1 , 1 ] 


The 2n— lifting of Lozano-Leal, supra, is equivalent to the specific 
choice of periodic windows of horizon length N=2n as shown in Eq. 
(21), which is a special case also to be excluded from the present 


invention . 
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N 



Py-\_0 , 0 f m . m f 0 f 1^1, ...fl] 


( 21 ) 


10 


15 


20 


Since each choice of p„ and Py gives rise to a unique lifting, there 
are a total of 2”^^ possible generalized liftings over a horizon of 
length N. Aside from the liftings of Eqs. (20) and (21), it appears 
that none of the other possible generalized liftings have been in- 
vestigated in the literature. It will now be shown that many of the 
other possible generalized liftings have very useful properties, 
particularly those of N>2n-1 . 

Zero Annihilation 

Focus will now be on those liftings for which the transmission 
zeros of the (squared down) lifted system lie at the origin. For 
notational convenience, we define the "small" matrix if, by, 

H^=SyHSj < 22 ) 


The quantity appears in many expressions and will play an 
important role in subsequent proofs. The matrix H, can be obtained 
directly from p„ and Py by writing p„ along the top of and Py along 
the side of H. This arrangement is depicted graphically in FIG. 3. 


25 
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The matrix is then the submatrix defined by the elements of H 
having I's along both borders. Two properties of interest 
concerning are. 

Output Tracking (OT) Condition: 

^ (23) 

Input Tracking (IT) Condition: 

hX-I (24) 

10 where I is the identity matrix and the superscript t denotes the 
Moore-Penrose inverse. It is noted that both the OT and IT 
conditions are satisfied if is square and invertible. The OT 
condition is intended for gain adjustments in control systems, and 
the IT condition is intended for equalization of communication 
15 channels. 

The property of placing transmission zeros of the lifted plant 
to the origin is characterized in the following; 

Lemma 1 {Zero Annihilation) Assume that windows p„ and py satis fy^ 
Zero Annihilation (ZA) Conditions: 


BSj=0 

(25) 

A(Sy)‘^=0 

(26) 


where f 


25 


S,=P7(pJ; Sy=W{Py)} s;=w{p^y) 



F94128 


20 


Then, 

(i) the generalized lifting of Eq. (19) has a simplified 
represen tat ion. 


, and 


(27) 


10 


15 


20 


Ys (k) SyASyY^(k-l) +SyHSjUg(k) (28) 

as illustrated in FIG. 4, and 

(ii) Y° is unobservabls from Y, and has stable (deadbeat) dynam- 
ics. 

Furthermore , 

(tat) at H, tn Eq. (23) is square and invertible, then the 
transmission zeros of the lifted transfer function in Eq. (19) from 
a. to Y. are annihilated (i.e., lie at the origin) . and 
(iv) it the OT condition Eq. (24) is satisfied then the trans- 
mission zeros of the "squared down" lifted transfer function Eq. 
(19) from V (where U.^h'.v) to Y, are annihilated. 

Proof: Results (1) and (li) follow by substituting the 

ZA conditions Egs. (25) and (26) into Eq. (19), to give. 




Ys(k-l) 

Ys'ik) 


Ys{k-1) 

Usik) 


Vsik-1) 




( 29 ) 
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where , 




Ys(}^)=C\ 


Ys{k-1) 

Usik-1) 





SyAS^ 0 0 


SyHSj 

^P= 

SyASy 0 0 

II 

SyHSj 


0 0 0 


I 


(30) 


(31) 


Cp=[s^s; 0 Oj; Dp=H^ (32) 

According to standard definitions (cf., E. J. Davison and S. H. 
Wang, "Properties and calculation of transmission zeros of linear 
multivariable systems," Automatica, Vol. 10, pp. 643-658, 1974) 
values of X satisfying. 


det 


|Ap-AJ Bp 


(33) 


are the transmission zeros of the transfer function from Ug(k) to 

Y,(k) defined by the state-space model (Ap,Bp,CpDp) . Consider the 
following identity, 

det ( G ) =det ( G 22 ) det ( - 0,^022 ) (34 ) 


where , 
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G& 


^ 11^12 

^ 21^22 


and the inverse of G,, is assumed to exist. Assuming that H, is 
invertible, the identity Eq. (34) can be applied to Eg. (33), which 
gives upon substituting Eqs. (31) and (32), 


-X.I 


det 


L 


X 

X 


0 0 
-X.I 0 


= (-X) 




0 -X.I 


which proves (iii) • Result (iv) follows by an identical analysis 
assuming that the OT condition Eq. (23) holds, and that the simpli- 
fied plant (AprBp,Cp,Dp) in Eqs. (29) -(32) has been squared down by 

a precompensator ffj . 

Lemma 1 is important since it gives conditions which and Py 
must satisfy for the generalized lifting G(p„, Py) to have its 
transmission zeros at the origin. All results in Lemma 1 can be 
simply understood by comparing the block diagram of FIG. 3 with 
FIG. 4 and noting all of the blocks that have vanished under the ZA 
conditions. It is seen that y/ no longer couples into the 
subsystem. Furthermore, the y/ subsystem has become deadbeat i.e. , 
all of the poles of the y/ subsystem are at the origin. Most 
importantly, there is now only one forward path from C7, to y,. 
Clearly something drastic has happened to the system zeros. Rigor- 
ously, if H 3 is square invertible, result (iii) of Lemma 1 states 
that the transmission zeros of the transfer function from U^(k) to 
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Ys(k) have been placed to the origin (i.e., annihilated). If is 
not square but the OT condition holds, result (iv) of Lemma 1 
states that the zeros the lifted plant "squared down" by a precoia- 
pensator Hg are annihilated. 

Extended Horizon Liftings 

The result of using extended horizons for the liftings 
introduces a new class of liftings which satisfy the conditions of 
Lemma 1 . 

Theorem 1 A class of generalized liftings G(p„, p^) which satisfy 
both the ZA and OT conditions is of the following form, 

Extended Horizon Lifting (OT Former 



Pu If •••flf p /1/1,...,1, 0,..a,0] 



q -1 n 



-A 



^ \ 
0 , . . . , 0 , 



(36) 


where elements m^O and £^0 are arbitrary, q=0, p e is an arbi- 
trary (or null) 0-1 vector chosen identically in both p„ and Py/ and 
n>0 is the order of the irreducible plant of the form given by Eq. 
(1). Furthermore, if the plant (system) of Eq, (1) is obtained by 
a zero-order hold (ZOH) digitization of a continuous- time system, 
the integer q^O can be chosen arbitrarily. 
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Proof-. By definition, the matrix B in Eq. (5c) has the sparse 
form B=[0JXJ where is a matrix of all ”0” elements, and 

By construction of S„ from p„ in Eq. (35), the nonzero 
elements of multiply only elements of in the product BS^ ■ 
5 Hence Bsj=0. Likewise, the matrix A in Eq. (5a) has the sparse 
form A = [OjXJ where is a matrix of all "0” elements, and 

By construction of Sy from Py (i.e., 0-1 complement of Py 
in Eq. (36)), the nonzero elements of Sy multiply only elements of 
O. in the product A(S/)^. Hence A(S/)^=0 and the ZA conditions of 

10 Eqs. (25) and (26) are satisfied. 

In order to show the OT condition, first consider the case 
where gsO. Then choice of Eqs. (35) and (36) ensures that the small 
matrix B. is of the form shown in FIG. 3, 


15 


i p n 

p(x\ Sr o\ 

® n U 1 -X’ 


(37) 


where ^ e R^"^ is lower triangular with a nonzero diagonal (and 
hence is invertible), and Ke B""" is given by. 


20 



■^2n+<l-l 



(38) 


Since ^is full rank, it follows from the special structure of Eq. 
(37) that Hg is full rank if IK has full rank. To show that IK has 
full rank, let (A,b,c) be any minimal (i.e., controllable and ob- 


25 
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servable) state-space realization of the transfer function Eg. (1) . 
The Markov parameters {h^} can be written as hj^=cA^~^b,i=l, . . . ,oo . 
Substituting into Eg. (38) gives. 




... cA^Jb 


... cA""^-^b 


( 39 ) 


3<=0a ‘^Cf 


( 40 ) 


where J is a reversed identity (i.e., J=[e„ . . . , e„] , f=[e,, . . . , e J ) , 
and O and C are observability and controllability matrices of 
(A,b,c), respectively. Since {A,b,c) is controllable and observ- 
able, it follows that 0 and C are each full rank. The fact the 
system or plant of Eg. (1) is obtained by a ZOH digitization im- 
plies that AeR"’"' is full rank (i.e., A is a state-transition 
matrix). These facts together imply that Min Eg. (37) is full 
rank and hence is invertible. 

If the system of Eg. (1) is not obtained by ZOH digitization, 
the matrix A may not be full rank. However, the results still hold 
wxth the restriction that g=0 since in this case is replaced by 
identity matrix I in relation Eg. (40), which is always full rank. 
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Theorem 1 is important because it extends Lozano's 2n-lifting 
to a much larger class of liftings which enjoy the same zero anni- 
hilation properties, and thus provides greater latitude (extra de- 
grees of freedom) in designing control systems. Note that the 
extended horizon liftings have total horizon length N=m+£+p+q+2n-l , 
which can be readily chosen longer than Lozano's lifting for which 
N-2n. Hence the phrase "extended liftings horizon". It will be seen 
that the extra degrees of freedom provided by N=k+2n-l, where 
k=m+£+p+q, for example, will overcome several difficulties 
associated with the 2n-lifting. For convenience, this new class of 
liftings is referred to herein as having N^2n—1 . 

The following result is essentially a "dual" to the previous 
theorem. 

Theorem 2 A class of generalized liftings G(p,,p^) which satisfy 
both the ZA and IT conditions is of the following form. 

Extended Horizon Lifting (IT Form) : 



N 

.A- 



0,...,0, p ,1,1, ...,1, o,...,o. 


n-1 

( * 

0 , ..., 0 ] 


( 41 ) 


m 

— » 

Py - [0,...,0, 





n-1 n 

^ ^ ■ 
Of mmmfOf 


( 42 ) 


where m^O and &>0 are arbitrary, q=0, peRF is an arbitrary (or null) 
0-1 vector chosen identically in both and p^; and n>0 is the 
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order of the irreducible plant Eg. (1). Furthermore, if the plant 
or system Eg. (1) is obtained by a zero-order hold (ZOH) 
digitization of a continuous-time system, the integer g^O can be 
chosen arbitrarily. 

Proof: The ZA conditions follow using the same arguments found 

in the proof of Theorem 1. Consider the IT condition Eg. (24). 
Given the lifting Eg. (41), has the form. 


Hs 


p n 

^ 0 ^ 

X H 

X X) 


(43) 


From the structure of Eg. (43) it follows that H, has full rank if 
both and IK have full rank. The remainder of the proof is identi- 
cal to the proof of Theorem 1. 


Note: Properties of Lozano's 2n-lifting (W=2n) follows directly 

from its interpretation as a special case of the liftings m2n-l in 

Theorem 1 and Theorem 2 . 

corollary 1 Lozano's 2 n-lifting Eg. (21) satisfies the ZA, OT and 
IT conditions. 

Proof: The lifting Eg. (21) is eguivalent to the special case of 

Theorems 1 and 2 where m=0, £=0 , p—lf p~[0]f g— 0 . 


the liftings in 
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ADVANTAGES 

The advantages of the extended horizon liftings in Eqs. (35) 

and (36), and in Eqs. (41) and (42) over Lozano's 2n-lifting are as 
follows : 

(i) If one chooses £>0 in Eqs. (35) and (36), there are more 
control inputs than outputs in the lifted system (i.e., a, > a^) . 
It is shown below with respect to application to plant inverse con- 
trol that these extra degrees of freedom can be used to design a 
controller which minimizes a quadratic control cost while simul- 
taneously satisfying a deadbeat tracking objective. This signifi- 
cantly reduces control gains compared with Lozano's lifting. 

(ii) If one chooses £>0 in Eqs. (41) and (42), there are more out- 
puts than inputs in the lifted system (i.e., a, > aj. it is shown 
below with respect to channel equalization in applications to com- 
munications that these extra degrees of freedom can be used to min- 
imize a quadratic error when estimating the input from measurements 
of the output. This is significant for reducing noise in problems 
of nonminimum phase channel equalization. Furthermore, the p vec- 
tor can contain additional message information to increase the 
channel throughput. 

(iii) If one chooses m>0 in Eqs. (35) and (36) or Eq. (37), there 
IS an extra m*T seconds of free time which can be used to perform 
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computations (where T is the sampling interval). Since m can be 
chosen arbitrarily, the use of extended horizon liftings for either 
control or equalization applications is not constrained by realtime 
computer limitations. This is particularly useful for adaptive 
implementations which involve additional computation. 

APPLICATION TO PLANT INVERSE CONTROL 

The placement of the transmission zeros to the origin by the 
class of extended horizon liftings of the form given by Eqs. (35) 
and (36) allows stable invertibility of the transfer function from 
Ug(k) to Yg(k) . A control law which will be discussed next deadbeats 
the response YJk) to follow the desired Y^(k) , subject to the mini- 
inization of a quadiratic control cost. 

To derive the desired controller, define the output error as, 

E(k)=Y^{k)-YJk) (44) 

Substituting Eq. (27) into Eq. (44) gives, 

E{k) =-SyASy^Y^(k-l ) -H^UJk) ^Y^(k) (45) 

Consider the problem of forcing the error in Eq. (38) to zero in a 

single step, while minimizing a quadratic control cost penalty, 
i.e. ^ 

min u:ik)U,(k) 


subject to 


F94128 


30 


E(k) =0. 


( 47 ) 


10 


In light of the output tracking condition of Eq. (23), thxs 
minimization problem can be solved uniquely (see D.S. Bayard and D. 
Boussalis, "Noncolocated structural vibration suppression using 
zero Annihilation Periodic control," 2nd IEEE Conference on Control 
Applications, Vancouver Canada, September 13-16, 1993) to give: 
Zero Annihilation Periodic (ZAP) Control Law: 


u°{k) =Hs (SyASyY^ ( k - 1 ) + ( k ) ) 


-K°YAk-l)^L "Y^(k) 


(48) 


(49) 


15 


where the corresponding feedback gains are defined as, 

K°^-H^SyASy 


(50) 




(51) 


Here the superscript "o" is chosen to emphasize the fact 
that the control nulls (i.e., deadbeats) the output. Also, in 
light of the OT condition, H. has full row rank and one can write 
(cf., Barnett, supra). 

For convenience the ZAP control law is summarized in the block 
diagram of FIG. 5. It provides the following result. 
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Lemma 2 (ZAP Control) Consider the closed-loop system arising 
from the OT extended horizon lifting Eqs. (35) and (36) under ZAP 
control Eg. (49 ) . Then, 

(i) The quadratic control cost Eq. (46) is minimized at each 
5 stage, subject to the deadbeat tracking constraint Eq. (47), 

(ii) All closed-loop poles are at the origin (i.e., the closed- 

loop response is deadbeat), and hence YJk) converges to YJk) in 
a single step. 

(Hi) The closed-loop system is internally stable (e.g., re- 

10 mains bounded) . 

Proof: Result (i) follows from the OT condition in Eg. (23) and 

well known minimum-norm properties of the Moore-Penrose inverse 

(cf., Barnett, supra). Now form the closed-loop system from the 

sxmplified lifted plant Eqs. (27) and (28) under ZAP control Eg. 
15 (49), 


Ysik) 


Ys(k~l) 

Ys(k) 


Ys(k-l) 

Usik) 


Us(k-l) 




20 where. 



{I-HJls)SyASj 0 

0 


HJIs 

Aal = 

s;(i-hsJh^s^)as/ 0 

0 

; Bci = 

S/HSX 


-HsSyASy 0 

0 


Hs 


(52) 


( 53 ) 
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Results (ii) and (iii) follow by substituting the OT condition 
Eg. (23) into Eg. (53), and noting that the resulting closed-loop 
matrix is stable with all of its eigen values at the origin. 

Result (i) of Lemma 2 is important because it indicates that 
control gains associated with using extended horizon liftings of 
the class denoted herein as N^2n-1 will be significantly reduced 
compared to those from using Lozano's lifting {N=2n) . Result (iii) 
of Lemma 2 is important because it ensures that the complementary 
output Yg remains "well behaved" even though it is not being 
controlled directly. 

Remark 1 Instead of deadbeat control, a pole placement scheme 
can be obtained by modifying the deadbeat constraint Eg. (47) to 
become E(k)=aE(k-l) in which case the ZAP control becomes 
U°{k) =K‘^YJk-l ) +L %(/c) -aL °E{k-l ) . 

Example 3 As an example, a 12-state nonminimum phase transfer 
function is shown in FIG. 6 (pole-zero plot, sampling time T=.025 
sec), adopted from the ASTREX flexible structure model(cf. Bayard, 
1993, supra and A. Das, J.L. Berg, G.A. Norris, D.F. Cossey, T.J. 
Strange III, and W.T. Schlaedgel, "ASTREX - A Unigue test bed for 
CSI research," Proc. 29th IEEE Conference on Decision and Control, 
Honolulu, Hawaii, pp. 2018-2023, December 1990). The open-loop 
response to an initial condition is shown in FIG. 7. A simulation 
xs fxrst run using the deadbeat control Eg. (49) with the Lozano's 
2n-lifting, m=0, 2=0, p=l , p=[0], n=12, g=0, N=2n=24, and using 7^=0 
(x.e., a vibration damping objective). 


The response is shown in 
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FIG. 8 to reach -3^10* at the output and 500 at the input. As 
expected from the theory, the response is deadbeat after a single 
horizon. However, this control law is unusable since the units are 
in volts, and the allowable range is only ±10 Volts. It is 
5 emphasized that these extraordinarily large responses are typical 
of the 2n-lifting due to the fact that the plant is inverted on a 
horizon of length N'T=0.6 sec. 

In an attempt to get a practical response with the 2n-lifting, 
the pole placement control of Remark 1 is used. The choice a=0.5 is 
10 made to get approximately a 10-second decay time (any slower would 
be worse than the open-loop response) . The results are simulated 
but not shown here since it turns out that the responses are re- 
duced 50%, and are still unacceptably large by several orders of 
magnitude. The ZAP control using an extended horizon lifting of 
15 the class N^2n-1 with m=0 , £=40, p=0, q=0, n=12 , N=63 was tried 
next. The results are shown in FIG. 9 where it is seen that both 
the input and output are well within the allowable ranges. The 
deadbeat nature of the response is also noted, as the vibrations 
are damped instantaneously after the first horizon at time T=3sec. 
20 This demonstrates the advantage of using the extended horizon lift- 
ings of the present invention instead of Lozano's 2n-liftings. 
Application to Channel Equalization 

An important problem in communications is that of equalizing 
(i.e., compensating for undesired cunplitude-frequency response) a 
25 nonminimum phase channel. The usual problem is that the channel 
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cannot be inverted in a stable fashion. However, using the input 
tracking (IT) extended horizon lifting form of Eqs. (41) and (42), 
this problem can be overcome. 

Referring to FIG. 10, let U^(k) be the sequence of messages to 
be sent by a transmitter T, and assume that an IT extended horizon 
lifting form given by Eqs. (41) and (42) is used to transmit the 
data (i.e., the signal sent is given by U{k) =sJUg{k) ) . The channel 
C is assumed to be a stable linear nonminimum phase transfer func- 
tion for a plant of the form given by Eq. (1) with a plant order n. 
Then from Theorem 2 and Lemma 1, the plant (channel) dynamics given 
by Eqs. (27) and (28) are 

Y^{k)=A^Y^{kl)^H^U^{k) ( 54 ) 

where 

A^=s^s;. ( 55 ) 

At the receiver R, the quantity Yg(k)=SyY{k) is measured, and it is 
desired to estimate the messages Ug(k) which were sent. For this 
purpose, an output prediction Yg{k) is formed as, 

Y,(k)=A^Y^{k-l)^H^UJk) ( 56 ) 

and an estimate is found by minimizing the least squares cri- 
teria. 


min(y^-yj^(r^-Y^) 


( 57 ) 
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Since the IT condition Eq. (24) holds for the extended horizon 
lifting, the unique solution to Eq. (56) is given by, 

Us = -HUA,YJk-l)-YJk)) (58) 

where one can write = since H, has full column rank 

(cf., Barnett, supra). It is noted that this dynamical system 
represents a stable inversion of the nonminimum phase channel 
characteristics. Furthermore, the use of an extended horizon 
lifting with £>0 has allowed channel inversion by least squares 
which provides smoothing in case of noise. Lozano's 2n-lifting in 
this application would not allow smoothing and would be very 
susceptible to noise. The use of p in the IT extended horizon 
lifting also provides a means for sending additional information 
with each packet, to improve the overall throughput. 

CONCLUSIONS 

A general class of liftings N>n has been denominated as N^2n-1 
in order to distinguish from the special case of Lozano (N=2n) and 
have been shown to have the same desirable zero annihilation 
properties of Lozano's 2n-lifting. In contrast to Lozano's 
2n-lifting, all the new liftings include horizons greater than 2n, 
i.e., they are of the extended horizon type N>2n-1 and o^>a^ or 
ay>a„. The use of extended horizons of the type N>2n-1 resolves 
many difficulties associated with Lozano's 2n-lifting. For 
example, a Zero Annihilation Periodic (ZAP) controller is defined 
for which the control gains can be significantly reduced as 
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coiRpaired to Lozano's lifting. This is due to a quadratic control 
cost which is minimized simultaneously with the deadbeat tracking 
objective. The effectiveness was shown in a simulation example 
where the control torque was reduced 4 orders of magnitude. 

5 As a dual result, it was shown that a related class of lift- 

ings enables the on-line equalization of nonminimum phase channels 
in communication systems. This overcomes the standard bottleneck of 
inverting the channel in a stable fashion. Here, channel inversion 
is accomplished by least squares estimation which provides smooth- 
10 ing in the case of noise. It is worth noting that this channel 

equalization approach can be made adaptive by using standard recur- 
sive algorithms since the new liftings are linear-in-the-para— 
meters . 

It is expected that the present invention (extended horizon 
15 liftings) will be useful in many areas of modern control, neural 

control, fuzzy control, adaptive control, communications, adaptive 
filtering, signal processing, or other applications where a stable 
system inverse is desired but not possible due to nonminimum phase 
constraints. 
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abstract 


periodic gain adjustment in plants of irreducible order, n, 
or for equalisation of communications channels is effected in 
such a way that the plant (system, appears to be minimum phase by 
choosing a horizon time N>n of liftings in periodic input and 
output windows and respectively, where i. is an integer 

chosen to define the extent (length) of each of the windows 
and p„ and n is the order of an irreducible input/output plant. 
The plant may be an electrical, mechanical or chemical system, xn 
which case output traolting (OT) is carrxed out lor feedback 
control or a communication channel, in which case input tracking 
(IT) is carried out. Conditions for OT are distinct from IT in 
terms of zero annihilation, namely lot OT and H, ^ 

IT, where the OT conditions are intended lor gaxn adjustments 
the control system, and IT conditions are intended for equalxza- 
tion for communication channels. 
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